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Abstract 

In this paper the massive neutrinos model arising due to the Snyder 
noncommutative geometry, proposed recently by the author is partially 
developed. By straightforward calculation it is shown that the masses 
of the chiral left- and right-handed Weyl fields treated as parameters 
fixed by experiments, lead to the consistent physical picture of the 
noncommutative geometry, and consequently yield renormalization of 
an energy of a relativistic particle and exact integrability within the 
proposed model. This feature of the model in itself both defines and 
emphasizes its significance and possible usefulness for both theory as 
well as phenomenology for high energy physics and astrophysics. 
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1 Introduction 

It was established recently by the author that the Dirac equation modified 
by the 7^-term arising due to the Snyder noncommutative geometry model, 
yields the conventional Dirac theory with nonhermitian mass, or equivalently 
to the massive neutrinos model given by the Weyl equation with a diagonal 
and hermitian mass matrix. The model describes 4 massive chiral fields re- 
lated to any original, i.e. non-modified, massive or massless quantum state. 
Due to spontaneous global chiral symmetry breaking mechanism it leads to 
the isospin-symmetric eff'ective field theory, that is composed chiral conden- 
sate of massive neutrinos. All these results violate CP symmetry manifestly, 
however, their possible physical application can be considered in a diverse 
way. On the one hand the effective theory is beyond the Standard Model, 
yet can be considered as its part due to the noncommutative geometry model 
contribution. On the other in the massive neutrinos model masses of the two 
left- and two right-handed chiral Weyl fields arise due to mass and energy 
of an original state, and a minimal scale {e.g. the Planck scale), and its 
quantum mechanical face becomes the mystic riddle. 

This paper is mostly concentrated on the quantum mechanics aspect. It 
is shown that the model in itself yields consistent physical explanation of the 
Snyder noncommutative geometry model and consequently leads to energy 
renormalization of an original quantum relativistic particle. Computations 
arising directly from the Schrodinger formulation of both the Dirac and the 
Weyl equations, are presented. First, the manifestly non hermitian modified 
Dirac Hamiltonian is discussed. Its integrability is formulated by straight- 
forward application of the Zassenhaus formula for exponentialization of sum 
of two noncommuting operators. It is shown, however, that this approach 
does not lead to well-defined solutions; for this case the exponents are still 
sums of two noncommuting operators, so that this procedure has a cyclic 
problem which can not be finished, and by this reason is not algorithm. For 
solving the problem instead of the Dirac equation we employ the Weyl equa- 
tion with pure hermitian mass matrix rewritten in the Schrodinger form. Its 
integration is straightforward and elementary. We present computations in 
both the Dirac and the Weyl representations of the Dirac gamma matrices. 

The paper is organized as follows. The Section 2 presents the motivation 
for further studies - the massive neutrinos model is recalled briefiy. Next, in 
the Sections 3 particle's energy renormalization is discussed. In the Section 

4 we present integrability problem for the modified Dirac equation. Section 

5 is devoted for the massive Weyl equation integrability, and the Section 6 
discusses some special case related to ultra-high energy physics. Finally in 
the Section 7 the results of the entire paper are summarized briefly. 
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2 The massive neutrinos model 

Let us recall briefly the massive neutrinos model resulting from [1]. The 
starting point is the noncommutative geometry model [2] of phase-space and 
space of a relativistic particle due to a fundamental scale i proposed by 
Snyder [3] (Cf. also Ref. [4]), and given by the following lattice model [5] 

x = ndx , dx = i , neZ — > i= -e^'"^ , lim £ = 0, (1) 

n n^oo 

where /q > is a constant, together with the De Broglie formula relating the 
coordinate x with its conjugate momentum p 

p=-. 2 

X 

Application of the Kontsevich star-product [6] to the phase space (a:,p) and 
two space points x and y 

x-kp = PX + Y^^(^^^ Cn{x,p), (3) 

x-ky = + Cn{x,y), (4) 



where for correctness a ~ 1, /3 are dimensionless constants, and Cn{f,g) are 
the Hochschild cochains, is yielding to the deformed Lie brackets (For review 
of deformation quantization see e.g. Ref. [7]) 

[x,pI = [x,p] + ^^(^^^ Bn{x,p), (5) 
[x,yl = [a:,?/] + ^(y) Bnix.y), (6) 

n=l ^ ^ 



where Bn{f,g) = Cn{f,g) — Cn{g,f) are the Chevalley cochains. Using 
[x,p] = —ih and [x,y] = 0, and doing the flrst approximation one obtains 

[x,p]^ = -ih+^Bi{x,p) , [x,y]^ = YBiix,y). (7) 

or in the Dirac "method of classical analogy" form [8] 

^[p,xl = l-^Bi{x,p) , ^[x,yl = ^Bi{x,y). (8) 
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Because, for any f,gE C°°{M) holds Bi{f,g) = 26{df A dg), one obtains 

^ [p, x]^ = l- ^{dx A dp) , ^ [x, y]^ = ^dx A dy, (9) 

where H in first relation was introduced for dimensional correctness. Applying 
now the lattice model (1) and the De Broglie relation (2) one receives 

l[^.Pi = i+^jy , ^[^'^L = -f^^ (10) 

that defines the Snyder model. Note that this model was studied in some 
aspect by previous authors [9], but the model is related to this direction to a 
slight degree. The model developed in this paper arise mostly from the idea 
of the papers [1]. 

If we consider £ as a minimal scale, e.g. Planck or Compton scale, then 
the model (10) can be rewritten in terms of the maximal energy e 

i|.,ri = i + lcV . = o(l) . (11) 

The lattice model (11) straightforwardly yield the contribution to the Ein- 
stein Hamiltonian constraint of Special Relativity 

_ cV ^ (7>^)2 = + ^cV , p, = [E,cp], (12) 

where m and E are mass and energy of a particle, and consequently leads to 
the generalized Sidharth 7^-term within the usual Dirac equation 

7^p^ ± mc^ ± icW) = , p^ = ih[do,cdi], (13) 

violating the Lorentz symmetry manifestly. In fact the equation (13) de- 
scribes 4 cases, that are dependent on the choice of the signs of mass m and 
the 7^-term. Here we will consider, however, positive signs case only. The 
negative ones are due to the changes e ^ — e and m —m in the results 
obtained from the positive signs case. 

Preservation of the Minkowski momentum space structure within the 
modified Einstein constraint (12) 

PmP^ = (7>m)' = (14) 

moves back considerations to the generic Einstein theory with e = 00, while 
application of the hyperbolic relation (14) within the modified Dirac equation 
(13) leads to two the conventional Dirac theories 

(7% + MiC^) ^± = 0, (15) 
4 
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where i/j"^ are the Dirac fields related to the nonhermitian mass matrices M±, 
that in general are dependent on an energy E and a mass m of an original 
quantum relativistic particle. It is of course the case of positive signs in the 
equation (13). In fact for any signs case there 2 the Dirac fields, so that the 
equation (13) describes 8 the Dirac fields. 

With using of the basis of projectors, M± can be decomposed as follows 

M± = f^%-^ + f^t^, (16) 
4 = (i ± \/'' - - 4^') ' (17) 



lit = -,[^± yje^ + 46mc2 - 4i?2) , (18) 

where /i^^ are the projected masses, or equivalently can be presented as a 
sum of its hermitian ^(M) and antihermitian 2l(M) parts 

M± =^(M±) + 2l(M±), (19) 
m±) = ^^^^U , 2t(M±) = ^^i^7^ (20) 

Introducing the chiral right- and left-handed Weyl fields ■0^ related to the 

1 ± 7^ 

Dirac field il)^ according to the standard transformation ^/'^^ = — ^ — 

one obtains two the massive Weyl equations with diagonal and hermitian 
mass matrices yU± 



, ^Ji 



± 



/i± 



± 







(21) 



so that we have totally 16 chiral fields describing by the Weyl equations (21) 
received from the Dirac equations (13). The massive Weyl theories (21) are 
the Euler-Lagrange equations of motion for the gauge field theory with chiral 
symmetry SU{'i)l^^ = SU{3)^ © SU{3)c 

C = C+ + C-, (22) 

where £^ are the Lagrangians associated with the fields i/j^mji^L as follows 

£± = i^h^p^^i^i + ijh^p,.iljt + fi^c^^Ri^^ + f^tc'^pL^t, (23) 

which is spontaneously broken to the composed gauge field theory with the 
isospin symmetry ^[/(2)™^ = SU{2)+ © SU{2)y 

C = V?"+(^Mp^ + ;,+ ^c2)^+ + f-(7% + ^;^/)V^-= (24) 
= (7% + ^e//c') ^, (25) 
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where fifjj are the effective mass matrices of the gauge fields ^/;^, and M^jf 



is the mass matrix of the effective composed field \i/ 



given by 



fi^^ff 



(26) 



The Lagrangian (25) describes the effective field theory - composed chiral 
condensate of massive neutrinos. 

In this paper we will consider both the Dirac equations (15) and the 
massive Weyl equations (21), but we are not going to discuss the gauge field 
theory (25) that will be studying in our next topical papers. We will assume 
that both the neutrinos masses (17) and (18) are real numbers, i.e. we will 
consider situation when the maximal energy e deforming Special Relativity 
is determined by a relativistic particle characteristics as follows 



-oo, — 2mc 



U 



U 



-2mc 



2mc^ 1 + \ 1 




U 



(27) 



or for the case of massless relativistic particle possessing an energy E is 

e G (-oo,-2|E|] U [2|^|,cx)). (28) 



3 Energy renormalization 

In fact existence of the massive neutrinos allows to explain in a consistent 
physical way the nature of the Snyder noncommutative geometry model. Let 
us see that by direct elementary algebraic manipulations the relations for 
masses of left- and right- chiral Weyl fields (17) and (18) can be rewritten in 
the form of the system of equations 

( fi^c^ + -V = - Aemc^ - AE^ 

2/ 2 (29) 
(^;,±c2-|) = e2 + 4emc2 - 4^2 

6 
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which allows to study dependence of the deformation energy parameter e 
and the particle energy E from the masses m, treated as physically 
measurable quantities. By subtraction of the second equation from the first 
one (29) one establishes the relation 



iemc 



(30) 



which after application of elementary algebraic manipulations allows to derive 
the deformation energy in the Snyder model (11) as 



1 - 



am 



(31) 



A maximal energy (31) does not vanish for all jjf ^ fxf ^ 0, and is finite for 
all /^R + /^L 7^ 8m. Here m is the mass of an original quantum state, and 
both fi^ and are assumed as physical quantities. In supposition all the 
masses can be fixed by experiments. In the case, when an original state is 
massless, one obtains 



e(m = 0) = (^if - ^if) = eo, 



(32) 



that is finite and non vanishing for finite /i^ 7^ and /x^ 7^ 0. In this manner 
we have 

'im 



1 + 



O 



5m 



for all l/i^ + /i^l > 8m, and 



8m 



O 



8m 



(33) 



(34) 



for all l/i^ + /i^l < 8m. On the other hand, however, addition of the second 
equation to the first one in (29) gives the relation 



2 (e^ - 4£;2) 



(35) 



which can be treated as the constraint for the energy of a relativistic 
particle, immediately solved with respect and presented in the canonical 
quadratic form with respect to the energy parameter e 



E^ 



e + 



1 ^ 






2 














3 J 





7 



(36) 
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By application of the deformation parameter energy (31) into the energetic 
constraint of a relativistic particle (36) one obtains the formula 



48 



8m \ 

8m 



7 + 



which for the case of originally massless state reduces into the form 



(37) 



3-4- 



f^tf^R 



{3t 



In fact, for given Eq the equation (38) can be used for establishment of the 
relation between masses of the neutrinos. In result one receives two possible 
solutions 

5 



± 4 i 

f^R = g/^L 



± Wl + 3 



En 



1^0 



(39) 



which are minimized for the value /ig = by the values 



(40) 



Interestingly, there is a possibility of the one solution between the masses 
and /i^ that is given by putting fiQ as a tachyonic mass 



and results in the relation 



(41) 



(42) 



For all l/i^ + /i^l > 8m the constraint (37) can be approximated by 

2" 



E^ - El 



± _ ± N 2 



8m 



o 



6m 



f^L+f^R 



(43) 



and for |/^^ + < 8m the leading approximation is 



E"^ -El 



± _ ± N 2 



5m 



4 2 /i± + /x^ 



O 



Qm 



, (44) 
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From the relation (35) one sees that, because the LHS as a sum of two squares 
of real numbers is always positive, it follows that the RHS must be positive 
also. In result we obtain the renormalization of a relativistic particle's energy 
E values 

-liEii, (45) 

Naturally, for the generic case of Special Relativity we have e = oo and by 
this energy E values are not limited. In this manner, in fact the Snyder 
noncommutative geometry results in energy renormalization of a relativistic 
particle. 



4 Integrability I: The Dirac equation 

The modified Dirac equation (15) can be straightforwardly rewritten in the 
Schrodinger evolutionary equation form (See e.g. the papers [11] and the 
books [12]) 

ihdoij^ = Hip^, (46) 
where in the present case the Hamilton operator H can be established as 

H = -zhn^d, - ^^ + ^^ cV + ^^"^^ cVt', (47) 

and splitted into its hermitian Sj{H) and antihermitian 2l(^) parts 

H = S){H)+^{H), (48) 

f,(i/) = -z^c/fS. - ^^^^cV, (49) 

Ql(H) = ^^"^^ cVt', (50) 
with (anti)hermiticity defined standardly 

j d^xtp^S){H)ij;^ = j d^xSj{H)tP^tP^, (51) 

d^xiJ^Ql{H)tp^ = - j d^x^{H)tp±tp^. (52) 

Note that in the case of equal masses fi^ = fxf = fi the antihermitian part 
(50) vanishes identically, so that the hermitian one (49) gives the full contri- 
bution, and consequently (48) becomes the usual Dirac Hamiltonian 

Hd = -7° {thcfd, + fic") . (53) 
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For this usual case, however, from (31) one concludes that 

e = 0, (54) 

so in fact the minimal scale becomes infinite formally i = oo, and by (37) 
relativistic particle's energy becomes E = i-^c^ with some mass /i. If we, 

however, take into account the tachyonic mass case fj. ^ ifj, = fx' then (53) 
becomes 

= [ihcfdi + ifi'c"^) , (55) 

and E = i^l^'^^- The relation (45), however, is not validate in this case. 

The full modified Hamiltonian (47) has nonhermitian character evidently, 
so consequently the time evolution (46) is manifestly non unitary. Its formal 
integration, however, can be carried out in the usual way with the following 
time evolution operator 



^±(x,t) = G(t,to)V^±(x,to) , G{tM)=Q^\^{-l i dTH{T)\ . (56) 



h 



to 



By this reason, the integrability problem for (46) is contained in the appro- 
priate Zassenhaus formula 

oo 

exp{^ + 5} = exp(A)exp(S) JJexpC„, (57) 

n=2 

C2 = -\C. (58) 

Cz = -^i2[C,B] + [C,A]), (59) 

C, = -^i[[C,A],A] + 3[[C,A],B] + 3[[C,B],B]), (60) 

where C = [A, B]. By the formula (56) one has identification 

A ^ A{t) = -^iyrS){H){r), (61) 



to 
t 



B = Bit) = -- / dT'ii{H){T), (62) 



to 

so that the commutator C is established as 

"f rt 



C=-^ [ dr' f dT"€{T',T"), (63) 

Jtn Jtn 



'to J to 

10 
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where 



9^{H){r')MH){r") . (64) 



Straightforward calculation of C can be done by elementary algebra 



2 

where we have applied the relations 

'J J ^ J = —J -J J -J , ^ -)/ 'y = —J J J J ^boj 

arising by employing the usual Clifford algebra of the Dirac matrices 

{^^,Y} = 2v'''U , {7',7^}=0 , f = z7V7V. (69) 

So, consequently one obtains the result 

itir',T") = 2^{H){r'MH)iT"), (70) 

that leads to the equivalent statement - for any times r' and r" the Poisson 
bracket between the hermitian Sj{H){t') and the antihermitian 2t(//)(r") 
parts of the total Hamiltonian (48) is trivial 

{m){T')MH){r")}=0. (71) 
Naturally, by simple factorization one obtains also 

C = 2AB , {A,B} = 0, (72) 



and consequently 



[C,A] = CA, (73) 

[C,B] = CB, (74) 

[[C,A],A] = 2[C,A]A (75) 

[[C,A],B] = 2[C,A]B, (76) 

[[C,B],A] = 2[C,B]A, (77) 
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and so on. In result the 4th order approximation of the Zassehnaus formula 
(57) in the present case is 

exp(y4) exp(i?) exp C2 exp C3 exp C4, 
1 



exp {A + B} 
C2 



C3 
C4 



iCA + 2CB), 



2 
1 

( CA' + 3CB' + 



(78) 
(79) 

(80) 
(81) 



For the case of constant in time masses and /i^ one determine the relations 



A 



C 



2h 

(/^L - f^jy 



{t - tof -ihcYdi + 



C 7 , 



(82) 

(83) 
(84) 



and consequently by elementary algebraic manipulations one establishes the 
Zassenhaus exponents as 



C2 
C3 



6^3 



X 



C4 



7' + (/^L - f^%y 



7 



(86) 



12^4 



X 



± _ ± N 2" 



7^ + 



(87) 



Exponents C„ show in a manifest way that the integrability problem is not 
well defined. Namely, the Zassenhaus coefficients C„ are still a sums of two 
noncommuting operators. The fundamental stage, i.e. the exponentialization 
procedure, must be applied again, so that consequently in the next step one 
has the same problem, i.e. the cyclic problem. Therefore this recurrence is 
not algorithm, that is the symptom of non integrability of (46). 
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5 Integrability II: The Weyl equation 



For solving the problem, let us consider the integrability procedure with 
respect to the massive Weyl equation (21). This equation can be straightfor- 
wardly rewritten in the form of the Schrodinger time evolution 



ihdo 



where the hermitian Hamilton operator H defines to the unitary evolution 











(89) 



so that the integration can be done in the usual quantum mechanical way. 
Integrability of (88) is well defined. The solutions are 



U{t,tu 



(90) 



where f/(t,to) is the unitary time-evolution operator, that for the constant 
masses is explicitly given by 



U{t,to) 



exp{--{t-to)H 



(91) 



and ■0^^(x,to) are the initial time to eigenstates with defined momenta 



(92) 



where the momenta and pf^ are related to the right- ijj^{x,to) or left- 
handed ip^{x,tQ) chiral fields, respectively. The eigenequation (92), however, 
can be straightforwardly integrated. The result can be presented in the 
compact form 



i^RA^^'to) = exp 



Pr,l {x-xo)ia 



h 



IpRA^O^to), 



(93) 



or after direct exponentialization 



I2 cos 



± 

Pr,l 



h 



[X - XQ)i 



— I 



■ ± 
P_RJ_^ 

h 



X - xo ha 



sm 



± 
Pr,l 

h 



X - Xo i 



± 

Pr,l 

h 



[X - XQ)i 



(94) 
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Currently the embarrassing problem presented in the integration procedure 
of the Dirac equation, discussed in the previous section, is absent. The 
Zessenhaus formula is not troublesome now because the matrix 7^ is by def- 
inition included into the Weyl fields, so that the Hamilton operator (89) is 
pure hermitian, and consequently the exponentialization (91) can be done 
is the usual way. At first glance, however, the mass matrix presence in the 
Hamilton operator (89) causes that one chooses at least two nonequivalent 
representations of the Dirac 7 matrices. Straightforward analogy to the mass- 
less Weyl equation says that the appropriate choice is the Weyl basis. On the 
other hand, however, the Hamilton operator (89) is usual hermitian Dirac 
Hamiltonian, so consequently the Dirac basis is the right representation. In 
this manner, in fact, we should consider rather both the chiral fields and the 
time evolution operator (91) strictly related the chosen representation (r) 

U'-it,to), (95) 
(^±J'-(x,to), (96) 
iij%JXxo,to) (97) 

where the upper index r = D,W means that the quantities are taken in the 
Dirac or the Weyl basis. The eigenequation (92), however, is independent 
on the representation choice, so that it physical condition - the fields have 
measurable momenta ^ . For full correctness, let us test both choices. 



U{t,to) 



5.1 The Dirac basis 

The Dirac basis of the gamma matrices is defined as 



7 



/ 
-/ 



Y 



a' 
-a' 



7 



I 

1 



(98) 



where / is the 2x2 unit matrix, and a* = [ax, (Jy, cTz] is a vector of the 2x2 
Pauli matrices 



1 

1 



-i 

1 



1 
-1 



(99) 



Consequently, by using of (98) the Hamilton operator (89) becomes 



H 



± 



(100) 
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and for the case of constant in time neutrinos masses yields a solution by the 
unitary time evolution operator U 



U 



D 



exp { -^-^(i - ^o) 



(101) 



After straightforward exponentialization (101) can be written in the compact 
form 



U 



D 



I 

/ 



COS 



c 



— I 



c 2 



X 



sm 



X 



4 + , f,^^.g^ 



(102) 



where we understand the all the functions are treated by the appropriate 
Taylor series expansions. 



5.2 The Weyl basis 

Equivalently, however, one can consider employing of the Weyl representation 
of the Dirac 7 matrices. This basis is defined as 



7 



' 


/ " 




' a' ' 




' -/ " 


/ 





-a' 




/ 



(103) 



For the choice of a representation in the form (10-3) the massive Weyl equation 
(88) is governed by the Hamilton operator (89) having the following form 



H 



(104) 



15 



L. A. Glinka / Energy renormalization and integrability within tlie massive neutrinos model 



Consequently, for the case of constant in time neutrinos masses one estab- 
lishes the unitary time evolution operator in the following formal form 



U 



w 



exp 



if, 

i—aOi 

c 



-i—aOi 
c 



(105) 



which after straightforward elementary exponentialization procedure can be 
presented in the form 



U 



w 



cos 



— I 



I 
/ 



h 



X 



sm 



X 



^ ~ ^0 ^2 / ,,±,,± 



h 



(106) 



Evidently, time evolution operator derived in the Weyl representation 
(106) has simpler form then its Dirac's equivalent (102). In this way the 
choices are not physically equivalent, i.e. will yield different solutions of 
the same equation. It is not, however, the strangest result. Namely, both 
the choices can be related to physics in different energy regions. So that it 
is useful to solve the massive Weyl equation in both mentioned representa- 
tions. It must be emphasized that strictly speaking the results obtained in 
this subsection are related to the massive Weyl equations presented in the 
Schrodinger time-evolution form (88). 



5.3 The space-time evolution 

Presently, one can employ the results received above, i.e. the momentum 
eigenequations (92), the spatial evolutions (93), and the unitary time evolu- 
tion operators (102) and (106), for an exact determination of the appropriate 
solutions of the massive Weyl equation (88) in both the Dirac and the Weyl 
representations of the Dirac gamma matrices. 
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5.3.1 Dirac-like solutions 

Applying first the Dirac representation, by elementary algebraic manipula- 
tions one receives straightforwardly the right-handed chiral Weyl fields in the 
following form 



cos 



sm 



±0^ 



sm 



■ ±0 



t-t 



h 



X exp < —I 



2h 



exp <j -^P%\x - xo)ia' \ (V'^)^ 



exp <! -^Pf{x - xo)^(^' \ {^i)o } X 



(t -to)}, 



(107) 



where for shorten notation {iP%i)q = ('0^^)'^(xo, to), = and 



\ 



±0' 



(108) 



Similarly, the left-handed chiral Weyl fields also can be determined in an 
exact way, the result is analogical to the right-handed case 



cos 



h 



+ 



sm 



+ ifi^c^- 



— (Pl ) 



sm 



■ ±0 



E^ipf) 

-^e^{pI ) 



X exp < —I 



E^{pf) 



exp <j - Xo)ia' \ {ijj^)^ 



exp <; -^P%\x - xo)^cT' \ (^|)^ }> x 



2h 



{t - to) 



(109) 
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5.3.2 Weyl-like solutions 

Similar line of thought can be carried out in the Weyl basis. An elementary 
calculation leads to the right-hand chiral Weyl fields in the form 



cos 



h 



sm 



■ ±0 

Wr c- 



t-t 



h 



±0\ 



sm 



exp <j -■ - Xo)ia' \ (ij^)^ + 



ew\-lpfi^-^oU'\i^t)^^ (110) 



where similarly as in the Dirac-like case we have introduced the shorten 
notation (V'^^^)?' = (^^^i)^(xo, to), /i± = /i± = f^ff^R and 



For the left-hand chiral Weyl fields one obtains the formula 



\ 



fill) 



LJ^E'^ipf) 



sm 



+ ipj^ c- 



-^E'^ip^ ) 



±0n 



sm 



'-^E-^ipf) 



E^ipf) 



exp <; -^pf^ix - Xo)ia' \ (V'l)^ + 



exp<i-^P^°(x-xo).an(^±)o^. (112) 



In this manner one sees that the difference between obtained solutions 
is crucial. Direct comparing of the Weyl-like solutions (110) and (112) with 
the Dirac-like solutions (107) and (109) shows that in the Dirac basis case 
there are different coefficients of cosinuses and sinuses, there is additional 
time-exponent, and moreover the functions M^{p^) and M^(p^°) having 
a basic status for both the solutions also have different form with respect to 
choice of the Dirac 7 matrices representation. 
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5.4 Probability density. Normalization 



If we know the chiral Weyl fields, then in the Dirac representation, one can 
derive the usual Dirac fields by the following way 



±\D 



12 
12 



12 
12 



where for shorten notation (^Z' ) — j , v j , anyj. ^ 
Similarly in the Weyl representation, the Dirac fields can be determined as 



(113) 



02 



02 

±\W- 

r) - 



(114) 



where also we have used the shorten notation {^J ) = (V' ) (x,t), and 
{''^^rlY^ — i^l^^RhY^ {^1 It is evident now, that in general these two cases are 
different from physical, mathematical, and computational points of view. In 
this manner, if we consider the quantum mechanical probability density and 
its normalization, we are forced to relate the Lorentz invariant probability 
density to the chosen representation 



fi^'^ = (f±)^'^(v;±; 



D,W 



d^xfi^''^ = I4. 
Using of (113) by elementary calculation one obtains 



D 



(115) 



(116) 



J-2 -L2 

2 2 

2 2 

(117) 



By application of (114) the probability density (115) becomes 



n 



w 



Oo 



02 



(118) 



Employing the normalization condition (116) in the Dirac representation one 
obtains the system of equations 



rf^a;(V7±^)^(V^±)^ + / d'xi^P^^n^t) 



(119) 
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which leads to 

= 1, 

In the case of Weyl representation one receives straightforwardly 



(120) 

(121) 
(122) 

(123) 
(124) 



In this manner one sees that in fact both the conditions (121), (122) and 
(123), (124)) are invariant with respect to choice of gamma matrices repre- 
sentations 

(125) 



D,Wrj.± \D,W 



that means they are physical. Using of the fact that full space-time evolution 
is determined as 



ii^tLf'^^M = t/^'^(t,to)(V^^,J^''^(x,to), 

one finds easily the condition 

^'^(^"^iJ,L)^'^(^,to)(^ij^'^(x,to) = 1. 



(126) 
(127) 

(128) 



By using of the spatial evolution (94) one obtains the relation 



\x — Xo\ 



-xsa sm 



± 

2-^fx - XQ)i 



h 



(129) 



where C = |('0^^)-^''^(xo,to)|^ is a constant, and — — 

nary part of the vector a\ The decomposition cTj = [ax, 0, az] + ^[0, —iay, 0] 
yields Qa'' = [0, —iay,0], and the equation (129) becomes 



it 



1l 



IS a imagi- 



C \ (rx\\ 



2 - 



|x — XqI 



-Gy sm 



± 
Prl 
2— ^(x - xo)i 



1. 



(130) 
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Introducing the change of variables (x — Xo)i x'i in the following way 

± 

= 2-y-(x - xo)i, (131) 



and the effective volume V due to the vector x' 



V% = I d'x' { U - tayx'y^-^^ \ , (132) 



\x' 



the equation (130) can be rewritten in the form 



= ^l2, (133) 



so that one obtains easily 

± 0\ 3/2 



{i^%rr''^{xo.t,) = I 2^ I -^e^v^e±, (134) 



where 9± are arbitrary constant phases. The volume (132) differs from the 
standard one by the presence of the extra axial (y) volume Vy 



sin I x' ' 




Vy = -iay J d^x'x'y ^"^^'^^ ' , (135) 

which is the axial effect and has nontrivial feature, namely 

on finite symmetrical spaces 
Vy = oo on infinite symmetrical spaces . (136) 
on sections of symmetrical spaces 

One sees now that the normalization is strictly speaking dependent on the 
choice of a region of integrability. For infinite symmetric spatial regions this 
procedure is not well defined, because the axial effect is infinite. However, one 
can consider some reasonable cases that consider the quantum theory on finite 
symmetric spatial regions. Moreover, the problem of integrability is defined 
with respect to choice of the initial momentum of the Weyl chiral fields 
In fact there are many possible nonequivalent physical situations connected 
with a concrete choice of this eigenvalue. The one of this type situations 
related to a finite symmetric spatial region, we are going to discuss in the 
next section as the example of the massive neutrinos model, which in general 
was solved in this paper. 
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6 The reasonable case 



Let us consider finally the reasonable case, that is based on the normalization 
in a finite symmetrical box and putting by hands the value of initial momenta 
of the chiral Weyl fields according to the Special Relativity 



± _ ± 



(137) 



For that simplified case the normalization discussed in the previous section 
leads to the following initial data condition 



in. 



± \D,W, 



Xq, to) 



where V = J cPx. Introducing the function E^{x,y) 



(138) 



E^ix,y)^c' 



X + y 



x^ 



(139) 



the right- and the left-hand chiral Weyl fields in the Dirac representation 
take the following form 



{i^^f{x,t) 



cos 



t - 1 



h 



1 



sm 



sm 



t-t 



h 



IC 



,±l 



exp -—fiji[x - xojiO- 



h 



ic 



exp -^IJ-l[x - Xo)i(r 



X 



exp < t 



0± ^ (^-^o)] 



(140) 



and 



{i;t)^{x,t) 



cos 



t - 1 



h 



+ 



t-t 



'E^if^lt.^) 



h 



IC 



exp <; - ^o)i(^ 



— I 



1 



sm 



^ ~ ^0 TT.D/,,± ,.±\ 



h 



IC 



exp <j -"-^jd%{x - Xo)ia' 



X 



exp < I 



0^ ^ (^-^o)] 



(141) 
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Similarly, introducing the function E {x,y) 



(142) 



for the studied case the right- and the left-hand chiral Weyl fields in the Weyl 
representation have a form 



and 



■ ± 2 



exp z6'-| 



sm 



■ ± 2 

+ \ — xrsm 



V 



t-t 



cos 



t-t 



h 



j^W( ,± ,,±^ 



t-t 



j^W( ,± 



h 



IC 



exp <j - Xo)icr 



exp 



ic 



fi^{x - xo)ia 




, (143) 



■ ± 2 



exp id J, 



sm 



■ ± 2 



V 



t-t 



cos 



h 



pW/,,± ,,±\ 



pW/,,± ,,±^ 



h 



exp 



zc 



/i^(x - Xo)jCr* } + 



IC 



exp <; -—iJji[x - xo)i(T 



■. (144) 



The "reasonable case" considered above is only the example following 
from the massive neutrinos model given by the massive Weyl equations (21) 
obtained due to the Snyder model of noncommutative geometry (11), and 
naturally it is not the only case. Actually there are many other possibilities 
for determination of the relation between the initial eigenmomentum values 
p^ ]^ and the masses /i^^ of the right- and left- hand chiral Weyl fields ip^^j^- 
However, the concrete choice (137) tested in this section presents a crucial 
reasonability contained in its special-relativistic-like character. Naturally 
this choice is connected with the special equivalence principle applied to 
the massive neutrinos in at the beginning of their space-time evolution, i.e. 
E 



± 

R,L 



P^l'c. This case, however, is also nontrivial from the 
high energy physics point of view [13], namely it is related to the region of 
ultra-high energies, widely considered in the modern astrophysics (See e.g. 
[14] and suitable references therein). So, the presented reasonable case of the 
massive neutrinos evolution in fact describes their physics in this region, and 
has possible natural application in ultra-high energy astrophysics. 
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7 Outlook 

In this paper we have discussed in some detail the consequences of the mas- 
sive neutrinos model arising due to the Snyder model of noncommutative 
geometry. The massive neutrinos model is a consequence of the Dirac equa- 
tions for a usual relativistic quantum state supplemented by the generalized 
7^-term [1]. In fact, Sidharth has suggested that this term could give a neu- 
trino mass, however, in spite of a good physical intuition he has finished 
considerations on a laconic statement only, with no any concrete calculations 
and propositions for a generation mechanism of neutrinos masses [15]. 

First we have considered the physical status of the Snyder model. By 
detailed calculation we have shown that, in contrast to Special Relativity 
theory, within the massive neutrinos model an energy of any original rela- 
tivistic massive or massless quantum state is strictly renormalized due to a 
maximal energy, directly related to a minimal scale £ being the deformation 
parameter in the Snyder noncommutative geometry. In this manner the Sny- 
der model has received a deep physical sense, that is in some partial relation 
to the Markov-Kadyshevsky approach [9]. 

Next the integrability problem of the massive neutrinos model was de- 
tailed discussed. First we have considered the modified Dirac equations, 
which rewritten in the Schrodinger form have yielded manifestly nonhermi- 
tian Hamiltonian being a sum of a hermitian and a antihermitian parts. By 
employing the 4th order approximation of the Zassenhaus formula we have 
proven that the procedure is not algorithm by the presence of the cyclic 
problem in exponentialization. Consequently, the Dirac equations are not 
integrable exactly. By this formal reason we have redefined the integrabil- 
ity problem with respect to the massive Weyl equations corresponding to 
the Dirac equations. The massive Weyl equations was also rewritten in the 
Schrodinger form, and by using of both the Dirac and the Weyl represen- 
tations of the Dirac 7 matrices, we have constructed its analytical exact 
solutions. We have shown that the normalizability of a solution is correctly 
defined only for special regions of spatial integration, i.e. sections of sym- 
metric spaces or finite symmetric spaces. 

Finally, the case related to ultra-high energy physics and astrophysics was 
shortly discussed. The all obtained results in general present interesting new 
physical content. In itself the investigated quantum-mechanical approach to 
the massive Weyl equation is novel. There are still possible applications of the 
proposed massive neutrinos model to phenomenology of particle physics and 
astrophysics, especially in the ultra-high energy region. The open question 
is the gauge field theory related to the massive neutrinos model, possessing 
some features of QCD. 
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